		
	* binscatter plot
	
		quietly: reghdfe math_std nothreat_cycle3 i.year i.sch_id $prevvars $gradevar $yeardosevar $experiencevar $sdemovars $yeardosevar  if subject == "Math", absorb(id) cluster(id#year) 
		local eq = `"Slope = "'
		local eq = `"`eq' `: display %5.4f _b[nothreat_cycle3]'"'
		local eq = `"`eq' (`:  display %5.4f _se[nothreat_cycle3]')"'
		binscatter math_std nothreat_cycle3 if subject == "Math", nquantiles(35) noaddmean control(i.year i.sch_id $prevvars $gradevar $yeardosevar $experiencevar $sdemovars $yeardosevar) absorb(id) ///
			mcolors(black) lcolors(gray) xtitle("Unmonitored Days (Residualized)") ytitle("Standardized Math Scores (Residualized)") b2title(`"`eq'"') xscale(r(-6 8)) xlabel(#8) ///
			savegraph("$output_dir\Figures\binscatter_math.png") replace
		
		
		quietly: reghdfe read_std nothreat_cycle3 i.year i.sch_id $prevvars $gradevar $yeardosevar $experiencevar $sdemovars $yeardosevar  if subject == "ELA", absorb(id) cluster(id#year) 
		local eq = `"Slope = "'
		local eq = `"`eq' `: display %5.4f _b[nothreat_cycle3]'"'
		local eq = `"`eq' (`:  display %5.4f _se[nothreat_cycle3]')"'
		
		binscatter read_std nothreat_cycle3 if subject == "ELA", nquantiles(35) noaddmean control(i.year i.sch_id $prevvars $gradevar $yeardosevar $experiencevar $sdemovars $yeardosevar) absorb(id) ///
			mcolors(black) lcolors(gray) xtitle("Unmonitored Days (Residualized)") ytitle("Standardized Reading Scores (Residualized)") b2title(`"`eq'"') xscale(r(-6 8)) xlabel(#8) ///
			savegraph("$output_dir\Figures\binscatter_reading.png") replace
		
		
				
			
		
		
		